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Abstract 


A  novel  time-frequency  signal  decomposition  approach  is  evaluated  for  the  detec¬ 
tion  of  maneuvering  air  targets  in  sea-clutter  using  High-Frequency  Surface-Wave 
Radar  (HFSWR).  This  approach  is  based  upon  time-frequency  analysis  and  its  ei¬ 
genvalue  decomposition.  Results  demonstrate  that  the  new  approach  provides  an  ef¬ 
fective  way  for  detecting  and  analyzing  maneuvering  air  targets  in  heavily  cluttered 
regions.  In  addition,  the  proposed  approach  successfully  detects  the  maneuvering 
target  in  all  cases  of  the  experiment.  In  contrast,  the  traditional  Fourier  transform 
method  fails  to  detect  air  targets  when  a  certain  degree  of  acceleration  is  present. 
This  is  because  of  the  phenomenon  known  as  Doppler  smearing.  It  is  statistic¬ 
ally  shown  that  the  new  approach  produces  satisfactory  results  even  with  very  low 
signal-to-noise  ratios.  The  results  obtained  are  robust  with  respect  to  the  algorithm 
parameters.  The  basis  of  the  method,  theoretical  derivations,  and  comparisons  with 
traditional  procedures  are  outlined  here-in. 


Resume 


Une  nouvelle  approche  de  la  decomposition  des  signaux  temps/frequence  est  eva- 
luee  en  vue  de  la  detection  de  cibles  aeriennes  en  manIJuvre  dans  un  clutter  de 
mer  au  moyen  d’un  radar  haute  frequence  a  ondes  de  surface  (HFSWR).  Cette  ap¬ 
proche  est  fondee  sur  V analyse  temps/frequence  et  la  decomposition  de  sa  valeur 
propre.  Les  resultats  montrent  que  la  nouvelle  approche  foumit  un  moyen  efficace 
de  detecter  et  d’analyser  les  cibles  aeriennes  en  manIJuvre  dans  des  zones  presen- 
tant  un  clutter  eleve.  En  outre,  l’approche  proposee  permet  de  detecter  la  cible  en 
manIJuvre  dans  tous  les  cas  de  l’experience.  Par  contre,  la  methode  classique  des 
transformations  de  Fourier  ne  permet  pas  de  detecter  les  cibles  aeriennes  lorsqu’un 
certain  degre  d’ acceleration  est  present.  Cela  s’explique  par  le  phenomene  qu’on 
appelle  trainage  Doppler.  On  montre  statistiquement  que  la  nouvelle  approche 
donne  des  resultats  satisfaisants,  meme  en  presence  de  rapports  signal/bruit  tres 
faibles.  Les  resultats  obtenus  sont  solides  en  ce  qui  conceme  les  parametres  des 
algorithmes.  Le  fondement  de  la  methode,  les  derivations  theoriques  et  les  compa- 
raisons  avec  les  procedures  classiques  sont  decrits  aux  presentes. 
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Executive  summary 


A  High-Frequency  Surface- Wave  Radar  (HFSWR)  is  an  effective  and  relatively 
low-cost  means  of  providing  over-the-horizon  surveillance  of  surface  vessels  and 
low-flying  aircraft  in  coastal  regions.  These  radar  systems  have  demonstrated  the 
ability  to  detect  and  track  surface  vessels  beyond  400  km  range  and  small  low-flying 
aircraft  out  to  120  km  range.  Thus,  these  systems  can  be  used  to  monitor  activity 
within  the  full  range  of  the  Exclusive  Economic  Zone  (EEZ).  The  technology  of¬ 
fers  highly  autonomous  24-hour-per-day,  7-day-per-week  operation  with  one  of  the 
lowest  operating  cost  per  unit  area  of  all  radar  types. 

The  Canadian  Department  of  National  Defence  (DND)  has  engaged  in  research  and 
development  of  HFSWR  technology  more  than  twenty  years.  One  of  the  objectives 
of  this  program  is  to  demonstrate  the  potential  of  employing  HFSWR  technology 
as  a  means  of  providing  wide-area  coastal  surveillance  of  surface  vessels  and  low- 
flying  air  targets  in  mid-latitude  and  arctic  regions.  An  operational  HFSWR  must 
be  able  to  operate  in  a  congested  signal  spectrum  using  limited  signal  bandwidth, 
and  in  the  presence  of  sea  clutter,  ionospheric  interference  and  man-made  noise. 
Before  an  operational  HFSWR  system  can  be  implemented,  it  must  demonstrate 
the  capability  in  a  wide  range  of  signal  and  physical  environments. 

One  of  the  dominant  sources  of  clutter  that  degrades  the  detection  capability  of  sur¬ 
face  and  low-flying  targets  using  HFSWR  is  sea-clutter.  Conventionally,  targets  are 
detected  by  the  Fourier  transform  or  Doppler  processing  method.  However,  the  per¬ 
formance  of  the  Fourier  method  degrades  if  the  target  has  high  acceleration.  For  an 
accelerating  target  the  Doppler  spectrum  smears.  The  degree  of  smearing  becomes 
higher  with  increasing  coherent  processing  interval  (CPI)  for  a  given  acceleration 
or  when  the  acceleration  is  increased  for  given  number  of  pulses.  If  the  smearing  is 
too  much,  then  the  Fourier  method  will  fail  to  detect  the  target. 

This  report  introduces  a  new  approach  for  the  detection  of  maneuvering  air  targets 
in  sea-clutter  using  the  time-frequency  analysis  and  its  eigenvalue  decomposition. 
Results  demonstrate  that  the  new  approach  provides  an  effective  way  for  detecting 
and  analyzing  accelerating  and  decelerating  air  targets  in  heavily  cluttered  regions. 
In  addition,  the  proposed  approach  successfully  detects  the  maneuvering  target  in 
all  cases  of  the  presented  experiment.  In  contrast,  the  traditional  Fourier  transform 
method  fails  to  detect  air  targets  whenever  the  acceleration  is  increased  for  a  given 
number  of  pulses.  This  is  because  of  the  phenomenon  known  as  Doppler  smearing. 
It  is  statistically  shown  that  the  new  approach  produces  satisfactory  results  even 
with  very  low  signal-to-noise  ratios.  The  results  obtained  are  robust  with  respect  to 
the  algorithm  parameters.  The  basis  of  the  method  and  comparisons  with  traditional 
procedures  are  outlined  here-in.  The  method  presented  here  is  not  restricted  to 
this  particular  application,  but  it  can  also  be  applied  in  various  other  settings  of 
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nonstationary  signal  analysis  and  filtering. 

Although  this  new  approach  demonstrates  the  improved  detection  performance  of 
an  accelerated  air  target  using  HFSWR  in  sea-clutter  environment,  a  more  detailed 
analysis  will  be  performed  in  a  future  study  in  relation  to  other  clutter  environments 
such  as  ionosphere,  meteors,  and  lighting  environments.  Furthermore,  the  s-method 
based  detectors  will  be  developed  for  real-time  operational  scenarios  and  will  then 
be  compared  with  existing  Fourier  Transform  based  detectors.  Once  the  new  ap¬ 
proach  outperforms  the  currently  available  detector  methods,  numerically  efficient 
algorithms  will  be  developed  for  real-time  operations. 


T.  Thayaparan,  L.  J.  Stankovic  and  M.  Dakovic  ;  2005;  A  novel 
approach  for  the  detection  of  maneuvering  air  targets  in  sea-clutter 
using  High-Frequency  Surface-Wave  radar;  DRDC  Ottawa 
TR  2005-251 ;  Defence  R&D  Canada  -  Ottawa. 
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Sommaire 


Le  radar  haute  frequence  a  ondes  de  surface  (HFSWR)  est  un  moyen  efficace  et 
relativement  peu  couteux  d’assurer  une  surveillance  au-dela  de  l’horizon  de  navires 
de  surface  et  d’avions  volant  a  basse  altitude  dans  les  regions  coheres.  Ces  systemes 
de  radar  ont  demontre  la  capacite  de  detection  et  de  poursuite  de  navires  de  surface 
au-dela  de  400  km  et  d’avions  volant  a  basse  altitude  jusqu’a  environ  120  km.  Ces 
systemes  peuvent  done  servir  a  la  surveillance  d’activite  dans  la  zone  economique 
exclusive  (ZEE)  au  complet.  La  technologie  offre  un  fonctionnement  hautement 
autonome  24  heures  sur  24,  7  jours  sur  7,  les  frais  de  fonctionnement  par  unite  de 
surface  etant  parmi  les  plus  faibles  pour  tous  les  types  de  radar. 

Le  ministere  de  la  Defense  nationale  (MDN)  du  Canada  est  engage  dans  des  tra- 
vaux  de  recherche  et  developpement  sur  la  technique  HFSWR  depuis  quinze  ans. 
L’un  des  objectifs  de  ce  programme  est  de  demontrer  la  possibility  d’employer  la 
technique  du  radar  HFSWR  comme  moyen  d’assurer  une  surveillance  cohere  eten- 
due  des  navires  de  surface  et  des  cibles  aeriennes  volant  a  altitude  moyenne  et  dans 
les  regions  arctiques.  Un  radar  HFSWR  operationnel  doit  pouvoir  fonctionner  a 
des  frequences  encombrees  en  utilisant  une  largeur  de  bande  limitee  de  signal  et  en 
presence  de  clutter  de  mer,  de  brouillage  ionospherique  et  de  bruit  d’origine  artifi- 
cielle.  Avant  qu’un  systeme  de  radar  HFSWR  operationnel  puisse  etre  mis  en  place, 
il  doit  demontrer  sa  capacite  dans  toute  une  gamme  de  milieux  de  signalisation  et 
de  fonctionnement. 

L’une  des  sources  dominantes  de  clutter  qui  degrade  la  capacite  de  detection,  au 
moyen  d’un  radar  HFSWR,  des  cibles  de  surface  et  des  cibles  volant  a  basse  altitude 
est  le  clutter  de  mer.  Normalement,  les  cibles  sont  detectees  par  la  methode  des 
transformations  de  Fourier  ou  de  traitement  Doppler.  Le  rendement  de  la  methode 
des  transformations  de  Fourier  a  cependant  un  effet  de  degradation  si  la  cible  a  une 
forte  acceleration.  Pour  une  cible  en  acceleration,  les  frequences  Doppler  presentent 
une  trainee.  Le  degre  de  trainage  s’accentue  avec  V augmentation  de  l’intervalle  de 
traitement  coherent  (ITC)  pour  une  acceleration  donnee  ou  lorsque  1’ acceleration 
est  accrue  pour  un  nombre  donne  d’impulsions.  Si  le  trainage  est  hop  eleve,  la 
methode  des  transformations  de  Fourier  ne  permet  pas  de  detecter  la  cible. 

Le  present  rapport  presente  une  nouvelle  approche  pour  la  detection  des  cibles  ae¬ 
riennes  en  manIJuvre  dans  un  clutter  de  mer  au  moyen  de  F  analyse  temps/frequence 
et  la  decomposition  de  sa  valeur  propre.  Les  resultats  montrent  que  la  nouvelle  ap¬ 
proche  assure  un  moyen  efficace  de  detecter  et  d’ analyser  des  cibles  aeriennes  en 
acceleration  et  en  freinage  dans  des  regions  presentant  un  clutter  eleve.  En  outre, 
l’approche  proposee  permet  de  detecter  la  cible  en  manIJuvre  dans  tous  les  cas  de 
l’experience  presentee.  Par  contraste,  la  methode  classique  des  transformations  de 
Fourier  ne  permet  pas  de  detecter  les  cibles  aeriennes  des  que  1’ acceleration  est 
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accrue  pour  un  nombre  donne  d’impulsions.  Cela  s’explique  par  le  phenomene 
qu’on  appelle  trainage  Doppler.  On  montre  statistiquement  que  la  nouvelle  ap- 
proche  donne  des  resultats  satisfaisants,  meme  en  presence  de  rapports  signal/bruit 
tres  faibles.  Les  resultats  obtenus  sont  solides  en  ce  qui  conceme  les  parametres  des 
algorithmes.  Le  fondement  de  la  methode  et  les  comparaisons  avec  les  procedures 
classiques  sont  donnes  aux  presentes.  La  methode  decrite  aux  presentes  n’est  pas 
limitee  a  cette  application  particuliere,  elle  peut  egalement  etre  appliquee  a  diverses 
autres  methodes  d’analyse  et  de  filtrage  de  signaux  non  stationnaires. 

Meme  si  la  nouvelle  approche  demontre  un  meilleur  rendement  en  ce  qui  conceme 
la  detection  d’une  cible  aerienne  en  acceleration  au  moyen  d’un  radar  HFSWR  dans 
un  milieu  presentant  du  clutter  de  mer,  une  analyse  plus  detaillee  sera  effectuee  dans 
une  etude  ulterieure  par  rapport  a  d’autres  milieux  presentant  du  clutter,  comme 
l’ionosphere,  les  meteores  et  les  milieux  lumineux.  En  outre,  les  detecteurs  fondes 
sur  la  s-methode  seront  mis  au  point  pour  des  scenarios  operationnels  en  temps  reel, 
avant  d’etre  compares  aux  detecteurs  fondes  sur  la  methode  des  transformations  de 
Fourier.  Une  fois  que  la  nouvelle  approche  offrira  un  meilleur  rendement  que  les 
methodes  actuelles  de  detection,  des  algorithmes  efficaces  sur  le  plan  numerique 
seront  formules  pour  les  operations  en  temps  reel. 


T.  Thayaparan,  L.  J.  Stankovic  and  M.  Dakovic  ;  2005;  A  novel 
approach  for  the  detection  of  maneuvering  air  targets  in  sea-clutter 
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1  Introduction 


One  of  the  dominant  sources  of  clutter  that  degrades  the  detection  capability  of  sur¬ 
face  and  low-flying  targets  using  High-Frequency  Surface  Wave  Radar  (HFSWR) 
is  sea-clutter.  In  a  coherent  HFSWR  radar,  the  dominant  form  of  sea-clutter  is  due 
to  scattering  from  sea  waves  with  wavelengths  which  are  half  of  the  wavelength  of 
the  radar  carrier  frequency  and  which  are  travelling  radially  towards  and  away  from 
the  radar  [1],  This  form  of  scattering  is  called  Bragg  scattering  and  it  corresponds 

to  the  wave  velocity  v  =  \J~fz,  where  A  is  the  radar  carrier  wavelength  and  g  is 
the  acceleration  due  to  gravity.  Thus  when  the  received  signal  of  a  pulse  Doppler 
radar,  taken  at  a  particular  range,  is  viewed  in  the  Fourier  transformed  or  Doppler 
domain,  we  can  see  two  sharp  peaks  at  frequencies  f Bragg  =  ±2u/A  = 

The  Bragg  components  of  the  clutter  are  called  1st  order  clutter.  In  addition  to 
the  Bragg  components,  the  Doppler  spectrum  contains  a  continuum  called  the  2nd 
order  clutter. 

Signals  are  commonly  analyzed  in  either  the  time  or  frequency  domain.  However, 
some  signals  exhibit  significant  time  variations  of  the  frequency  content.  For  these 
cases  time-frequency  representations  can  be  used,  since  they  combine  time  and  fre¬ 
quency  domain  analyses  to  yield  a  more  revealing  picture  of  the  temporal  localiza¬ 
tion  of  signals  spectral  components  [5,  14,  23,  24],  The  oldest  and  the  most  widely 
used  time-frequency  representation  is  the  short-time  Fourier  transform  (STFT).  In 
order  to  improve  its  concentration,  various  quadratic  representations  have  been  in¬ 
troduced  [3,  6,  12],  The  most  prominent  member  of  this  class  of  representations  is 
the  Wigner  distribution  (WD).  Inversion  properties  of  the  Wigner  distribution  and 
synthesis  of  a  signal  from  a  given  time-frequency  representation  have  been  studied 
[4,8,9,11,13], 

The  time-frequency  representation  referred  to  as  the  S-method  has  a  property  that, 
under  certain  assumptions,  its  value  for  multicomponent  signals  is  equal  (or  close) 
to  the  sum  of  the  Wigner  distributions  of  individual  signal  components  [17],  It 
has  been  also  used  as  a  model  in  the  implementation  of  time-scale  representations, 
time -varying  spectra  estimation,  detection  and  realization  of  higher  order  repres¬ 
entations  [2,  10,  15,  16,  18].  In  this  report,  the  S-method  is  used  to  introduce  a  new 
method  for  the  decomposition  of  multicomponent  signals  by  using  eigenvalues  and 
eigenvectors  of  an  appropriately  formed  matrix. 

The  proposed  decomposition  method  is  applied  in  the  analysis  of  HFSWR  radar 
signals.  Conventionally,  targets  are  detected  by  the  Fourier  transform  or  Doppler 
processing  method.  However,  the  performance  of  the  Fourier  method  degrades  if 
the  target  has  high  acceleration.  For  an  accelerating  target  the  Doppler  spectrum 
smears.  The  degree  of  smearing  becomes  higher  when  the  number  of  pulses  is 
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increased  for  a  given  acceleration  or  when  the  acceleration  is  increased  for  given 
number  of  pulses  [25].  If  the  smearing  is  too  great,  then  the  Fourier  method  will  fail 
to  detect  the  target.  The  case  of  highly  accelerating  targets  corresponds  to  the  ana¬ 
lysis  of  signals  with  fast  time  variations  of  the  frequency  content.  Therefore,  these 
kind  of  signals  should  be  analyzed  by  the  time-frequency  representations  rather 
than  by  the  Fourier  transform.  The  time-frequency  based  decomposition  can  also 
provide  the  extraction  of  individual  signal  components.  This  approach  will  be  effi¬ 
cient  in  separating  the  target  signal  from  an  undesirable  sea-clutter.  This  property  is 
of  a  particular  significance  in  the  HFSWR  signals  where,  besides  the  target  signal, 
there  is  a  clutter  that  is  primarily  due  to  scattering  from  the  surface  waves. 

The  results  obtained  by  applying  the  time-frequency  analysis  based  decomposition 
method  show  that  one  can  improve  the  detection  performance  of  the  HFSWR,  as 
well  as  add  new  insights  into  the  interpretation  and  processing  of  radar  signals,  with 
respect  to  the  traditional  Fourier  transform  methods  currently  used  by  HFSWRs. 

The  paper  is  organized  as  follows.  After  an  introduction,  the  theoretical  develop¬ 
ment  of  the  S-method  from  basic  principles  is  given  in  Section  2.  The  inversion 
of  the  Wigner  distribution  is  presented  in  Section  3.  A  review  of  the  S-method 
definition  and  its  basic  properties  is  given  in  Section  4.  This  representation  is  used 
for  defining  a  new  decomposition  method  of  multicomponent  signals  in  Section  5. 
Next,  the  signal  decomposition  is  discussed  in  Section  6  from  the  point  of  view  of 
a  target  signal  in  a  strong  sea  clutter.  Experimental  data  analysis  and  results  are 
presented  in  Section  7,  proves  the  efficiency  of  the  proposed  method,  and  illustrates 
improvements  in  readability  and  detection  of  the  target  signal.  The  noise  analysis  in 
the  algorithm  for  signal  decomposition  is  given  in  Section  8.  Conclusions  are  given 
in  Section  9,  Finally  the  Appendix  provides  a  high  resolution  form  of  the  S-method. 
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2  S-method:  Theory 


The  Fourier  Transformation  (FT)  of  the  signal  x(t)  is  defined  by: 


OO 


— oo 


(1) 


If  we  have  the  FT  of  the  signal  x(t),  we  can  calculate  the  signal  in  time  by: 


x(t ) 


—  [  X{u)ejutdu 

— OO 


(2) 


For  a  signal  x(t)  of  the  form: 

x(t)  =  e~jUTt  (3) 

its  FT  is  X(oS)  =  2k5{l o  —  u>t )•  Both  representations  of  the  signal  x(t),  in  the 
time  and  frequency  domains,  are  depicted  in  Figure  1.  The  reason  as  to  why  FT  is 
not  always  good  tool  for  representing  signals  is  illustrated  in  Figure  2  where  two 
different  signals  in  time  are  shown  in  the  frequency  domain. 

Although  the  phase  characteristic  of  FT  contains  the  information  about  the  time 
instants  when  the  spectral  components  occur  it  is  very  difficult  to  use  this  inform¬ 
ation.  As  a  result  we  can  see  in  Figure  2  two  signals  that  are  different  in  time,  but 
have  the  same  amplitude  spectral  representations.  Hence,  new  representations  have 
been  developed  whose  main  goal  is  to  provide  the  time  distribution  of  the  spectral 
content.  The  research  area  that  deals  with  these  representations  and  their  properties 


<D 
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Figure  1:  Representations  of  the  signal  x(t)  in  the  time  and  frequency  domains 


is  the  time-frequency  signal  analysis. 
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Figure  2:  Representations  of  two  different  signals  in  the  time  (left  hand  side) 
frequency  domains  (right  hand  side) 
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2.1  Short-Time  Fourier  Transform  (STFT) 

The  most  frequently  used  time-frequency  representation  is  the  Short-Time  Fourier 
Transform  (STFT).  The  idea  behind  the  STFT  is  to  introduce  a  window  function 
w(t)  that  truncates  the  analyzed  signal  x(t)  and  then  calculates  the  FT  of  the  trun¬ 
cated  signal.  The  STFT  is  obtained  by  sliding  the  window  along  the  signal.  The 
mathematical  formulation  of  the  STFT  in  the  analog  form  is: 


STFT(t,  uj)  =  J  x(t  +  r)w (r)e-jujT dr  (4) 

—  OO 

where  x(t  +  r)  is  the  signal  shifted  in  time  and  w(r)  is  a  window  function.  The 
illustration  of  the  signal  localization  for  the  STFT  calculation  is  given  in  Figure  3. 

The  discrete  form  of  the  STFT  is  given: 


OO  OO 

STFT(n,k)  =  w(m)x(n  +  m)e~:i^mk  =  w{m)xn{m)e~^mk  (5) 

m=— oo  m=— oo 


An  energetic  version  of  the  STFT  is  called  spectrogram.  It  is  defined  by: 
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SPEC(t,  cj)  =  \STFT(t,  lu)  |2  (6) 

The  following  gives  examples  of  the  STFT  for  different  signals  with  different  win¬ 
dow  widths. 

2.1.1  Example  1. 

Consider  a  complex  sinusoidal  signal  x(t)  with  constant  frequency: 


x(t)  =  e~j207Tt  (7) 

The  representation  of  the  signal  x(t)  in  the  time  domain  is  given  in  Figure  4.  The 
STFT  of  this  signal,  at  the  time  instant  t  =  0,  for  various  window  widths  (denoted 
by  T)  is  given  in  Figure  5. 

The  spectrogram  is  centered  at  the  signal  frequency  u  =  —20ir  =  —62.8.  The 
spectrogram  width  is  determined  by  the  width  of  the  FT  of  the  window.  Since  the 
Hanning  window  is  used,  if  we  assume  that  the  width  of  its  FT  is  equal  to  the  width 
of  the  main  lobe,  we  get: 


Wwindow  =  Air/T  (8) 

For  wide  windows,  when  T  — >  oo  ,  the  spectrogram  tends  to  a  delta  function  at 
u o  =  —62.8,  as  shown  in  Figure  5. 
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Figure  5:  The  spectrogram  of  the  signal  x(t),  at  the  time  instant  t  =  0  and  for 
various  window  widths  T.  Note  that  N  is  the  number  of  samples  within  window. 
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Figure  6:  The  representation  of  a  LFM  signal  x(t)  in  the  time  domain 


2.1.2  Example  2. 

Consider  a  Linear  Frequency  Modulated  (LFM)  signal  x(t): 


x(t)  =  e~jl607Tt2  (9) 

The  representation  of  the  signal  x(t)  in  the  time  domain  is  given  in  Figure  6.  The 
STFT,  at  the  time  instant  t  —  0,  for  different  window  widths  T,  is  given  in  Figure  7. 
The  window  function  w(t)  used  in  this  example  and  all  other  examples  is  a  Hanning 
window  function. 

Unlike  the  previous  example  where  the  STFT  is  dependent  only  on  the  window 
width,  in  this  example  the  STFT  is  dependent  on  the  window  width  and  non¬ 
stationary  feature  of  the  signal.  To  explain  the  above  mentioned,  we  will  introduce 
the  concept  of  instantaneous  frequency.  The  first  derivative  of  the  phase  is  referred 
to  as  the  instantaneous  frequency: 


ujft)  =  <j/(t) 

The  instantaneous  frequency  of  the  considered  signal  x(t)  is: 


(10) 


Ujft)  =  — 3207T  t  (11) 

If  t  is  within  the  range  —  ^  <t<  j,  it  follows  that  coft)  will  be  in  the  range: 


320t iT  .  .  320v rT 

—  -  Wi(f)  £  — 


(12) 


The  width  of  the  STFT  of  the  signal  x(t)  is  dependent  on  the  width  of  the  window 
and  the  width  of  the  signal: 
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Figure  7:  The  spectrogram  of  the  Linear  Frequency  Modulated  Signal  x(t)  in  the 
time  instant  t  =  0  and  for  various  window  widths  T.  N  is  the  number  of  samples. 
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fftot  kb window  T  signal  (13) 

Since  we  are  using  the  Hanning  widow  as  a  window  function  w(t),  it  follows  that 
Wwindow  =  ^,  and  from  that  Wsignai  =  3207 tT,  resulting: 

Wtot=4^+  320ttT  (14) 

It  indicates  that  the  optimal  window  width  could  be  obtained  from  dWtot / dT  =  0 
as  Toptimai  =  1/8.94.  The  spectrogram  calculated  using  this  widow  width  is  close 
to  the  one  depicted  in  Figure  7  with  T  =  1/8  .  Note  that  the  optimal  window  width 
is  highly  signal  dependent.  When  the  analyzed  signal  is  unknown,  the  optimal 
window  width  also  unknown.  This  is  always  the  case  in  practice. 

The  spectrogram  width  for  a  small  T,  for  example  T  =  1/32,  is  determined  by  the 
first  term  Wwindow  =  47 t/T  =  1287T  =  402  (see  Figure  7).  For  a  wide  window,  for 
example  T  =  1,  the  spectrogram  width  is  determined  by  the  second  term  Wsignai  = 
3207 rT  =  3207T  =  1005  (see  Figure  7). 

In  order  to  show  that  the  same  conclusions  hold  for  any  considered  time  instant, 
in  Figure  8  we  have  presented  the  spectrogram  at  the  time  instant  t  =  —0.2.  The 
two-dimensional  spectrogram  for  all  considered  time  instants  —0.5  <  t  <  0.5  is 
shown  in  Figure  9. 


2.2  Wigner  Distribution  (WD) 

In  order  to  efficiently  process  the  linear  frequency  modulated  signal,  the  Wigner 
Distribution  (WD)  is  used.  The  mathematical  formulation  of  the  Wigner  distribu¬ 
tion  in  the  analog  form  is: 


WD(t,u)=  /  w{^)w(~)x(t  +  ^)e  JUJTdr  (15) 


In  the  discrete  domain  the  WD  reads: 


OO 

WD(n,k)=  w(m)w(— m)x(n  +  m)x*(n  —  m)e~^mk  (16) 

m=— oo 


10 


DRDC  Ottawa  TR  2005-251 


Figure  8:  Spectrogram  of  Linear  Frequency  Modulated  Signal  x(t)  at  time  instant 
t  =  —  0.2  and  for  various  window  widths  T,  N  is  the  number  of  samples. 
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Figure  9:  The  spectrogram  of  the  Linear  Frequency  Modulated  Signal  for  various 
window  widths  N  (number  of  samples) 
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2.2.1  Example  3. 

Consider  the  same  LFM  signal  x(t)  as  the  one  denoted  by  formula  (9). 

x(t)  = 

The  representation  of  the  signal  x(t)  in  the  time  domain  is  given  in  Figure  10.  The 
Wigner  distribution  of  this  signal  is  calculated  in  the  subsequent  manner: 


WD(t,u)  =  J  we(T)e-jl60^t+^2ejl60^t-^2e-juTdr=  (17) 

—  OO 

OO  OO 

=  J  we{r)e-j320wtTe-juJTdr 

—  OO  — OO 

=  We(u  +  3207rt) 

where  we{r )  =  w(^)w(— |)  and  We(oj )  =  FT{we(r)}. 

Thus  we  can  conclude  that  the  Wigner  distribution  of  the  linear  frequency  modu¬ 
lated  signal  behaves  in  the  same  way  as  the  spectrogram  of  a  signal  with  constant 
frequency.  The  position  of  the  Wigner  distribution  is  determined  by  the  instantan¬ 
eous  frequency,  while  its  width  is  equal  to  the  width  of  the  FT  of  the  window.  It 
tends  to  a  delta  function  for  very  wide  windows.  In  order  to  illustrate  this  conclu¬ 
sion  we  presented  the  Wigner  distribution  at  the  instant  t  =  0  for  various  window 
widths  in  Figure  11.  Comparing  Figure  7  and  Figure  11,  the  advantage  of  the 
Wigner  distribution  over  the  spectrogram  for  LFM  is  obvious. 


J  we(r)e-jiuj+3207Tt)TdT  = 
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Figure  11:  The  Wigner  Distribution  of  the  Linear  Frequnecy  modulated  Signal  at 
the  instant  t  =  0,  for  various  window  widths  T  and  the  number  of  samples  N  within 
the  window 


We  have  presented  the  same  results  at  another  instant  t  =  —0.2  in  Figure  12.  For 
all  instants  the  Wigner  distribution  in  the  time-frequency  domain  is  presented  in 
Figure  13. 


2.2.2  Example  4. 

Consider  two  LFM  signal  components  x(t): 

x{t)  =  Xl(t)  +  x2(t)  =  e +  e-jl2nft*+jl0nft 

The  representation  of  the  signal  x(t)  in  the  time  domain  is  given  in  Figure  14. 
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Figure  12:  The  Wigner  Distribution  of  the  Linear  Frequnecy  modulated  Signal  at 
the  instant  t  =  —0.2,  for  various  window  widths  T  and  the  number  of  samples  N 
within  the  window 
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Figure  13:  The  Wigner  Distribution  of  the  Linear  Frequency  Modulated  Signal  for 
various  window  widths  N  (number  of  samples) 


Figure  14:  The  representation  of  two  LFM  signal  components  x(t)  in  the  time 
domain 
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The  Wigner  distribution  of  this  signal  is  calculated  in  the  subsequent  manner: 


WD(t,u)  = 


we(r)  xi(t  +  +  x2(t  +  xi(t  -  +  x2(t  -  ^)1  e  ja,Tdr 


—  OO 
OO 


we(r)xi(t  +  -  ^)*e  JUTdT  + 


+  /  uy(r)x2(t  +  ^)x2(t-  ^)*e  JUJTdr  + 


+2Re{  /  we(r)a:i(t  + ^)x2(t- ^)*e  JWTdr} 


where  we(r)  =  w(^)w(— |)  and  FFe(u;)  =  TT{we(r)}. 

From  the  last  formula  we  see  that  besides  the  Wigner  distribution  of  auto-terms: 

OO  OO 

/  we(r)xi(t  +  7j)xi(t  —  ^)*e~juTdT  and  J  we(r)x2(t  +  \)x2(t  —  ^)*e~juTdr, 

—  OO  — OO 

OO 

there  exists  a  cross-term:  2Re{  J  we{r)xi{t  +  ^)x2(t  —  ^)* e~ju>T dr} . 


We  calculated  the  Wigner  distribution  of  two  signal  components  for  various  window 
widths  at  the  time  instant  t  =  0.  It  is  shown  in  Figure  15.  We  see  that  besides  two 
auto-  terms  there  exists  a  strong  cross-term. 


The  same  results,  at  another  instant  t  =  0.035,  are  given  in  Figure  16.  The  value 
of  the  cross-term  changes  over  time  due  to  its  oscillatory  nature  can  be  seen  from 
Figures  15  and  16.  For  all  instants  the  Wigner  distribution  in  the  time-frequency 
domain  is  presented  in  Figure  17. 


2.3  S-Method 

In  order  to  eliminate  the  cross-terms  a  class  of  so-called  reduced  interference  dis¬ 
tributions  is  proposed.  All  of  them  are  based  on  the  Cohen  definition  of  the  quad¬ 
ratic  time-frequency  distributions,  with  a  kernel  function  being  a  low  pass  two- 
dimensional  function.  These  distributions  are  in  fact  forms  of  two-dimensional 
smoothed  Wigner  distribution.  By  smoothing  the  Wigner  distribution,  oscillatory 
cross-terms  are  reduced.  In  reduced  interference  distributions,  two  important  prob¬ 
lems  arise  and  should  be  stressed: 

First:  By  smoothing  the  Wigner  distribution  the  auto-terms  are  also  smoothed  [22], 
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Figure  15:  The  Wigner  Distribution  of  two  LFM  signal  components  x(t)  at  the 
instant  t  =  0,  for  various  window  widths  T  and  the  number  of  samples  N  within  the 
window 
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Figure  16:  The  Wigner  Distribution  of  two  LFM  signal  components  x(t)  at  the 
instant  t  =  0  .035,  for  various  window  widths  T  and  the  number  of  samples  N  within 
the  window 
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Figure  17:  The  Wigner  Distribution  of  two  Linear  Frequency  Modulated  Signal 
components  for  various  window  widths  N  (number  of  samples) 
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Second:  The  calculation  complexity  is  very  high,  including  need  for  oversampling 
in  most  of  the  cases. 

Here,  we  will  present  a  computationally  very  simple  method,  that  can  produce  the 
same  auto-terms  as  in  the  Wigner  distribution,  but  without  cross-terms.  It  will  be 
referred  to  as  the  S-method. 

The  S-method  is  derived  from  the  relationship  between  the  STFT  and  the  Wigner 
distribution,  which  reads: 


WD(t,u) 


oo 

■  J  STFT(t,  LO  +  9)  ST  FT*  (t,  cu  -  9)d9 

—  OO 


(19) 


Proof:  The  Wigner  distribution  definition  reads: 


WD(t,  u)=  /  w(^)w(-^)x(t  +  ^)x*(t-^)e  JUJTdr 


From  the  STFT  definition 


we  get 


OO 

STFT(t,u)  =  J  x(t  +  r)w(T)e-juTdr 


x(t  +  t)w(t)  =  —  /  STFT(t,  9)ej9Tdd 


Replacing  this  into  the  Wigner  distribution  definition  we  get: 


WD(t,u)  =  J  w{-l)x*(t~l) 


T  J  STFT(t,  9)ej6T/2d9 


e~jojrdT 


=  &  I  STFT(t,  9) 


J  w{—  7;)x*(t  —  |)e  d/2',Tdr  d9 


with  —  |  =  A  we  get 


=  b  f  STFT (f.  0) 


—  OO 
OO 


J  w(X)x(t  +  X)e  9^2XdX  d9 


=  l  f  STFT(t,  9)  ST  FT*  (t,  2  u  -  9)d9 

— OO 

with  9  =  lu  +  9  follows  : 

OO 

=  X  f  STFT(t,u  +  9)STFT*(t,uj-9)d9 
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This  ends  the  proof. 

A  discrete  version  of  the  previous  relation  reads: 

N/2 

WD(n,  k)  =  STFT(n,k  +  i)STFT*(n,k-i)=  (20) 

i=-N/2 

STFT(n,  k  +  i)STFT*(n,  k  -  i)  1 


N/2 


\STFT(n,  k)\2  +  2Re  E 


2=1 


The  mathematical  formulation  of  the  S-Method  (SM)  in  the  discrete  form  is: 


L 

SM(n,  k)  =  J2  STFT(n,k  +  i)STFT*(n,k-i)  = 

i=—L 

=  | STFT{n,  k) |2  +  2R e{STFT(n,  k  +  l)STFT*(n,  k  -  1)}+ 
+2R  e{STFT(n,  k  +  2)  ST  FT*  (n,  k  -  2)}  +  ... 
+2Re{5TFT,(n,  k  +  L)STFT*(n,  k  -  L)} 

where  L  determines  the  number  of  terms  used  in  the  calculation. 

It  can  also  be  written  in  the  form: 

N/2 

SM(n,  k)  =  P(i)STFT(n,  k  +  i)STFT*  (n,  k  -  i) 

i=—N/2 

where  P(i)  =  1  for  i  \  <  L  and  P(i)  =  0  for  other  values  of  i. 

The  S-method  with  L  terms  can  be  written  in  the  form: 

L 

SML(n,  k)=Y^  STFT(n ,  k  +  i)STFT*(n,  k  -  i) 

i=—L 

In  particular,  for  given  values  of  L  we  have: 

For  L  =  0,  the  S-method  is  identical  to  the  spectrogram 

SM0(n,  k)  =  \STFT(n ,  k) |2  =  STFT(n,  k)STFT*(n ,  k) 


For  L  =  1 

SMi(n,  k)  =  \STFT(n,  k)\2  +  2Re{STFT(n,k  +  l)STFT*(n,k  -  1)} 
=  SM0(n,  k )  +  2Re{STFT(n,  k  +  1  )STFT*(n,  k  -  1)} 
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For  L  =  2 

SM2(n,k )  =  \STFT(n,  k)\2  +  2Re{STFT(n,  k  +  l)STFT*(n,  k  -  1)}+ 
+2Re{STFT(n,k  +  2)STFT*(n,k  +  2)} 

=  SMi(n,  k)  +  2Re{STFT(n,  k  +  2 )STFT*(n,  k  -  2)} 

and  so  on. 

The  illustration  of  the  S-method  vector  calculation  for  L  =  2  is  given  in  Figure  18. 

Note  that  SM  with  L  terms  is  obtained  by  adding  one  more  term  to  the  S-method 
with  L-l  terms: 

SML(n,  k)  =  SML_i(n,  k)  +  2Re{STFT(n,  k  +  L)STFT*(n,  k  -  L)} 

Comment  on  the  number  of  terms  L: 

The  S-method  will  produce  the  same  auto  terms  as  the  Wigner  distribution  if  we 
take  L  such  that  (2 L  +  1)  is  equal  to  the  auto  terms  width  in  discrete  domain  (i.e., 
to  the  number  of  samples  within  the  auto  term).  In  practice  it  means  a  few  terms, 
for  example  L  =  3, 4, ...,  8, 9,  since  most  of  the  auto  term  energy  is  located  around 
its  maximal  value.  We  will  see,  in  the  next  examples,  that  the  performance  of  the 
S-method  is  insensitive  to  L  values  in  a  wide  range  of  L  values.  A  precise  mathem¬ 
atical  proof  can  be  given  that  the  S-method  can  produce  the  Wigner  distribution  for 
each  component  separately,  in  those  regions  of  the  time-frequency  plane  where  the 
components  do  not  overlap. 

2.3.1  Example  5. 

Consider  the  same  LFM  signal  x(t)  denoted  by  formula  (9).  The  representation  of 
the  signal  x(t)  in  the  time  domain  is  given  in  Figure  19. 


x  (t)  =  e~i160wt2 

We  calculated  SM  of  a  LFM  x(t)  at  the  instant  t  =  0,  fixed  window  width  N  =  128, 
and  different  L  values.  Results  are  shown  in  Figures  20  and  21.  For  all  instants  the 
S-Method  in  the  time-frequency  domain  is  presented  in  Figure  22. 

2.3.2  Example  6. 

Consider  two  LFM  signal  components  x(t)  denoted  by  the  formula  (18). 

X(t)  =  XX(t)  +  X2(t)  =  +  e-jl2Kft*+jlOKft 
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Figure  18:  The  illustration  for  the  S-method  calculation  for  L  =  2:  Signal  (first  row), 
STFT  calculated  as  the  signal’s  fft  (second  row),  STFT  shifted  for  one  sample  left 
and  right  (third  and  fourth  rows),  STFT  shifted  for  two  samples  left  and  right  (fifth 
and  sixth  rows).  The  S-method  is  equal  to  the  sum  of  the  squared  second  row  + 
double  real  part  of  the  product  of  third  and  fourth  rows  +  double  real  part  of  the 
product  of  fifth  and  sixth  rows. 
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Figure  19:  The  representation  of  a  Linear  Frequency  Modulated  Signal  x(t)  in  the 
time  domain 


15000 


Figure  20:  SM  of  a  Linear  Frequency  Modulated  Signal  x(t)  at  the  instant  t  =  0, 
fixed  window  width  N  =  128,  and  different  L  values 
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Figure  21:  SM  of  a  Linear  Frequency  Modulated  Signal  x(t)  at  the  instant  t  =  0, 
fixed  window  width  N  =  128,  and  different  L  values 


26 


DRDC  Ottawa  TR  2005-251 


frequency  frequency 

Figure  22:  SM  of  a  Linear  Frequency  Modulated  Signal  x(t),  with  fixed  window 
width  N  =  128,  and  different  L  values 
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Figure  23:  The  representation  of  two  LFM  signal  components  x(t)  in  the  time 
domain 


The  representation  of  the  signal  x(t)  in  the  time  domain  is  given  in  Figure  23. 

We  calculated  SM  of  two  LFM  Signal  components  x(t)  at  the  instant  t  =  0,  fixed 
window  width  N  =  128,  and  different  L  values.  Results  are  shown  in  Figures  24 
and  25.  In  Figure  26  only  positive  values  for  SM,  for  L=0...12  values,  are  shown. 
The  S-Method  for  all  instants  in  the  time-frequency  domain  is  presented  in  Figure 
27. 

2.3.3  Example  7. 

Consider  three  LFM  signal  components  x(t)  denoted  by  formula: 


X(t)  =  Xx(t)  +  X2(t)  =  e-j8nft2~jl5nft  +  +  e~j8nft*+j15nft  (21) 


The  representation  of  the  signal  x(t)  in  the  time  domain  is  given  in  Figure  28. 

We  calculated  SM  of  three  LFM  signal  components  x(t)  at  the  instant  t  =  0,  fixed 
window  width,  N  =  128,  and  different  L  values.  The  result  is  shown  in  Figure 
29.  In  Figure  30  only  positive  values  for  SM,  for  L=0...12  values,  are  shown.  The 
S-Method  for  all  instants  in  the  time-frequency  domain  is  presented  in  Figure  31. 

Figure  32  shows  the  S  method  of  three  Linear  Frequency  Modulated  Signal  com¬ 
ponents  x(t)  with  different  signal  rates,  with  fixed  window  width  N  =  128,  and 
various  L  values.  Figure  33  shows  the  S  method  of  three  nonlinear  frequency  mod¬ 
ulated  signal  components  x(t),  with  fixed  window  width  N  =  128,  and  various  L 
values. 
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Figure  24:  SM  of  two  Linear  Frequency  Modulated  Signal  components  x(t),  with 
fixed  window  width  N  =  128,  and  different  L  values 
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Figure  25:  SM  of  a  two  component  Linear  Frequency  Modulated  Signal  x(t) ,  with 
fixed  window  width,  N  =  128  ,  and  different  L  values 
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Figure  26:  Only  non-negative  SM  of  two  Linear  Frequency  Modulated  Signal 
components  x(t),  with  fixed  window  width  N  =  128,  and  L=0...  12  values 
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frequency  frequency 


Figure  27:  SM  of  two  Linear  Frequency  Modulated  Signal  components  x{t),  with 
fixed  window  width  N  =  128,  and  different  L  values 


Figure  28:  The  representation  of  three  LFM  signal  components  x(t)  in  the  time 
domain 
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Figure  29:  SM  of  three  component  Linear  Frequency  Modulated  Signal 
components  x(t),  with  fixed  window  width  N  =  128,  and  different  L  values 
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Figure  30:  Only  non-negative  SM  of  three  Linear  Frequency  Modulated  Signal 
components  x(t),  with  fixed  window  width  N  =  128,  and  different  L  values 
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Figure  31:  SM  of  three  Linear  Frequency  Modulated  Signal  components  x(t),  with 
fixed  window  width  N  =  128,  and  different  L  values 
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Figure  32:  The  S  method  of  three  Linear  Frequency  Modulated  Signal  components 
x(t)  with  different  signal  rates,  with  fixed  window  width  N  =  128,  and  various  L 
values 
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Figure  33:  The  S  method  of  three  nonlinear  frequency  modulated  signal 
components  x(t),  with  fixed  window  width  N  =  128,  and  various  L  values 
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3  Wigner  distribution  based  inversion  and 
decomposition 


A  discrete  form  of  the  Wigner  distribution  is  defined  by 

N/2 

WD(n,  k)  =  f(n  +  m)f*(n-m)e-j^2mk,  (22) 

m=—N/2 

where  we  assume  that  the  signal  f(n)  is  time  limited  within  \n\  <  N/2  and  omit 
a  constant  multiplication  factor  of  2.  Inversion  relation  for  the  Wigner  distribution 
reads 

N/2 

f(n  +  m)f*(n  -  m)  =  ^  WD(n,  k)ej^im(2k) . 

k=-N/2 

After  substitutions  n\  =  n  +  m  and  n2  =  n  —  m  we  get 

N/2 

(23) 

k=—N/2 

For  cases  when  (rti  +  n2)/ 2  is  not  an  integer,  an  appropriate  interpolation  is  per¬ 
formed  in  order  to  calculate  WD{{ni  +  n2)/ 2,  k )  Introducing  the  notation, 

N/2 

R(nun2)  =  j^  ^  WD( (24) 

k=—N/2 


we  get 

R{ni,n2)  =  f{ni)f*(n2).  (25) 

Matrix  form  of  (25)  reads 

R  =  f(n)f*(n),  (26) 

where:  f  (n)  is  a  column  vector  whose  elements  are  the  signal  values,  f*(n)  is  a  row 
vector  (Hermitian  transpose  of  f  (n))  and  R  is  a  matrix  with  the  elements  R(ni,n2), 
defined  by  (24). 

'For  better  understanding  of  the  calculation  procedure  note  that  relation  (23)  is  a  discrete  coun¬ 
terpart  of  the  Wigner  distribution  inversion  in  analog  domain,  that  reads: 

1  r°° 

^  j  WD((t i  +  t2) /2,u)e°Ni-^) du. 

By  discretizing  angular  frequency  to  =  kAuj  and  time  t\  =  niAt,  t2  =  n2At.  with  appropriate 
definition  of  discrete  values,  assuming  At  =  1,  we  easily  obtain  (23). 
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As  for  any  square  matrix,  the  eigenvalue  decomposition  of  R  reads 

JV+l 

R  =  QAQt  =  ^2  A,Uj(n)u*(n),  (27) 

2=1 

where  A*  are  eigenvalues  and  u,(n)  are  eigenvectors  of  R.  By  comparing  (26)  and 
(27),  it  follows  that  the  matrix  with  elements  of  form  (24)  can  be  decomposed  by 
using  only  one  non-zero  eigenvalue.  Note  that  the  energy  of  eigenvectors  is  equal 
to  1,  by  definition  ||u1(n)||2  =  1.  By  comparing  (26)  and  (27),  having  in  mind 
that  there  is  only  one  non-zero  eigenvalue,  we  have  f(n)f*(n)  =AiUi(n)if[(n)  = 
(v^Ui(n))(>Aiui(ri))*  or 


Ai 


N/2 

iif(«)ii2=  E  f»  =  Ei . 

n=—N/2 


resulting  in 


A i  =  Ef8(i  -  1), 


(28) 


where  8(i)  denotes  Kronecker  symbol.  Eigenvector  ui(n)  is  equal  to  the  signal 
vector  f  (n)  up  to  the  constant  amplitude  and  phase  factor.  Therefore,  an  eigenvalue 
decomposition  of  the  matrix,  formed  according  to  (24),  can  be  used  to  check  if  an 
arbitrary  2D  function  D(n,  k )  is  a  valid  Wigner  distribution. 


The  same  relations  can  be  used  in  signal  synthesis.  We  start  from  a  given  func¬ 
tion  D(n,  k ),  calculate  (24)  and  perform  eigenvalue  decomposition  (27).  The  first 
(largest)  eigenvalue  and  corresponding  eigenvector  produce  a  signal  such  that  its 
Wigner  distribution  is  the  closest  possible  Wigner  distribution  to  the  given  arbitrary 
function  D(n,  k ),  [11]. 


Now,  this  property,  along  with  the  S-method,  will  be  used  for  signal  decomposition. 
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4  S-method  definition  and  basic  property 


Note  that  a  definition  of  the  STFT  is 

N/2 

STFT(n,  k)  =  f{n  +  m)e~j^ ^mk .  (29) 

m=—N/2 

Its  relationship  with  (22),  as  derived  in  [17],  is 

N/2 

WD(n,k)  =  STFT(n,k  +  l)STFT*(n,k-l). 

l=-N/2 

This  relation  has  led  to  the  S-method  definition  [17,  21]: 

L 

SM(n,  k)  =  ^Yl  STFT(n ,  k  +  l)STFT*{n,  k  -  l)  (30) 

l=—L 


SM(n,  k)  —  P(l)STFT(n,k  +  l)STFT*(n,k-l )  (31) 

l=—N/2 

with  the  window  function  P(l )  =  1/(N  +  1)  for  |/|  <  L  and  P(l)  =  0  elsewhere. 
The  S-method  can  produce  the  representation  of  a  multicomponent  signal  such  that 
the  distribution  of  each  component  is  its  Wigner  distribution,  avoiding  cross-terms. 


Proposition :  Consider  a  multicomponent  signal 

M 

f(n)  = 

i=i 


where  fi(n)  are  monocomponent  signals.  Assume  that  the  STFT  of  each  component 
lies  inside  the  region  Dfn,  k),  i  =  1,2, ...,  M.  Denote  the  length  of  7-th  region 
along  k,for  a  given  n,  by  277;  (n),  and  its  central frequency  by  kofn).  The  S-method 
of  f(n )  is  equal  to  the  sum  of  the  individual  Wigner  distributions,  W Dfn,k), 
i  =  1, 2, ...,  M,  of  each  signal’s  component, 

M 

SM(n,  k)  =  J2  WDfn,  k),  (32) 

2=1 


if  the  regions  Dfn,  k),  i  =  1,2, ...,  M,  do  not  overlap,  Dfn,  k )  fl  Dj(n,  k)  =  0  for 
i  f  j,  and  the  number  of  terms  L  in  (30),  for  a  point  (n,  k),  is  defined  by: 


L(n,  k ) 


F>i(n)  —  | k  —  k0i(n)\  for  (n,k)  G  Dfin^k) 

0  elsewhere . 


(33) 
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The  proof  of  the  proposition  is  very  similar  to  the  one  provided  for  the  continuous 
S-method  case  and  can  be  found  in  [19]. 

Note  1 :  Real  M- component  signals  may  be  considered  as  2M-component  complex 
signals  with  each  region  Di(n,  k )  being  associated  with  the  region  Di+M  (n,  —k). 
Thus,  there  is  no  need  for  removing  negative  frequency  components  in  real  signals 
[20], 

Note  2:  Any  window  in  (30)  with  a  constant  number  of  terms  L  >  m‘dxnk{L(n,  k )} 

M 

produces  SM(n ,  k)  —  ^2  WD^n,  k ),  if  the  regions  D^n,  k ),  *  =  1,2, M,  are  at 

1=1 

least  2 L  apart  along  the  frequency  axis,  i.e.,  | k0i(n)  —  k0j  (n) \  >  B,  (n)  +  Bt  (n)  + 
2 L,  for  each  i,  j  and  n. 

This  is  the  S-method  with  constant  value  of  L,  as  it  was  originally  introduced  in 
[17,  18],  The  signal  dependent  method  (33)  would  be  more  accurate,  but  also  more 
complex.  Constant  number  of  terms  L  is  used  here  in  numerical  realization  since 
it  is  much  simpler  for  implementation,  producing  satisfactory  and  robust  results. 
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5  Decomposition  of  multicomponent 
signals 

For  each  signal  component  fi(n)  we  can  write  its  inversion  formula,  corresponding 
to  (23),  as 

N/2 

Tm  E 

k=-N/2 

i  =  l,2,...,M, 

if  we  knew  the  Wigner  distribution  WD^n,  k)  of  this  component.  By  summing  the 
above  relations  for  i  =  1, 2, M  we  get 

M 

= 

2=1 

N/2  M 

T^i  E  E  WDi(r^p-,ky*'Hni~m)- 

k=-N/2  i= 1 

The  equation  (32),  for  the  signals  that  satisfying  the  presented  conditions,  this  rela¬ 
tion  reduces  to: 

M 

Y Mnl)fi(n2 )  = 

i= 1 

N/2 

TO  E  (34) 

k=-N/2 

By  denoting 

N/2 

RSM(n1,n2)  =  ^  Y  (35) 

k=—N/2 

and  using  the  eigenvalue  decomposition  of  matrix  with  the  elements 

RsM(ni,n2),  we  get 

iV+1 

RsiW  =  Y  ^iui(n)ui(n)- 
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As  in  the  case  of  Wigner  distribution,  we  can  conclude  that  A,  =  Eft ,  i  =  1 , 2 , . . . ,  M 
and  Aj  =  0  for  i  —  M  +  1, N,  i.e., 

M 

k  =  EflS(i-l).  (36) 

i=i 

The  eigenvectors  u j(n)  will  be  equal  to  the  signal  components  up  to  the 

phase  and  amplitude  constants,  since  the  components  orthogonality  is  assumed  by 
the  Proposition.  Amplitude  constants  are  again  contained  in  the  eigenvalues  Aj. 
Thus,  the  reconstructed  signal  can  be  written  as 

M 

frecfa)  ^  ^  (^0 

1=1 


It  is  equal  to  the  original  signal,  up  to  the  phase  constants  in  each  component.  When 
we  have  several  components  of  different  energies  f]  (n),  f2  (n), ....  fM(n)  and  when 
they  are  all  of  equal  importance  in  analysis,  we  can  use  normalized  values  of  the 
signal  components  and  calculate  the  time-frequency  representation  of 

M 

f nor(n)  =  y^fe(A)Uj(n) 

1=1 

by  using  the  weights  k{ A)  =  1  in  the  signal,  i.e.,  by  using  the  eigenvectors  as  signal 
components  (Capon’s  way  of  weighting  in  the  minimal  variance  high  resolution 
approach,  the  Appendix). 

When  there  exists  a  very  strong  disturbing  signal  such  as  sea-clutter  in  the  HFSWR 
signal  we  can  omit  the  first,  strongest  component,  and  define  the  reconstructed 
signal  as 

Mi 

frec{n)  = 

i= 2 

where  Mi  is  the  expected  number  of  components. 


5.1  Illustrative  example 

Consider  a  signal  whose  analog  form  reads: 

x(t )  = 

k= 2 
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Figure  34:  Spectrogram  and  its  eigenvalue  decomposition  (upper  row);  Wigner 
distribution  and  its  eigenvalue  decomposition  (bottom  row). 


within  the  interval  —128  <  t  <  127,  where  uj2  =  —  cv3  =  —  |,  oq  =  — 
6U5  =  f ,  cu6  —  §,  co7  =  d2  =  d7  =  0,  d3  =  d5  =  -64  and  d4  =  d6  =  64. 

The  sampling  interval  is  At  =  1.  The  spectrogram,  calculated  according  to  (29), 
is  presented  in  Fig.34,  upper  row.  The  Wigner  distribution  is  presented  in  Figure 
34,  bottom  row.  Based  on  the  Wigner  distribution,  the  elements  of  matrix  R  are 
calculated  by  using  (24).  Eigenvalue  decomposition  (27)  of  this  matrix  produces 
exactly  one  non-zero  eigenvalue,  Ai  =  390.92  (A2  =  0.00,  A3  =  0.00,  ....),  being 
equal  to  the  total  signal  energy  Ex  =  390.14  (within  the  numerical  calculation 
error),  as  expected  from  (24)-(28). 

To  illustrate  that  eigenvalue  decomposition  of  the  spectrogram  (that  is  the  first  step 
in  S-method  calculation)  does  not  produce  a  meaningful  result,  we  have  repeated 
this  procedure  by  using  the  spectrogram  instead  of  the  Wigner  distribution  in  cal¬ 
culation  of  the  elements  of  R,  according  to  (24).  The  eigenvalue  decomposition 
obtained  in  this  way  proves  that  the  spectrogram  can  not  be  related  to  a  sum  of  the 
Wigner  distributions  of  the  signal  components,  Figure  34,  upper  row  (right) . 

The  S-method  of  the  same  signal  is  calculated  by  using  (29)  and  (30)  with  L  = 
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12.  The  obtained  results  are  depicted  in  Figure  35.  The  matrix  R sm  is  formed 
according  to  (35).  Its  eigenvalue  decomposition  results  in  the  same  number  of 
non-zero  eigenvalues  as  the  number  of  signal  components.  Eigenvalues  correspond 
to  the  components  energies,  while  the  eigenvectors  correspond  to  the  normalized 
signal  components,  up  to  the  phase  constants.  The  time-frequency  representation 
of  the  eigenvectors  is  shown  in  Figure  36.  First  seven  components  correspond  to 
the  signal,  while  the  remaining  ones  are  with  very  small  eigenvalues.  Energies  of 
discrete  signal  components  are:  Ex  =  119.40,  E2  =  50.13,  E?>  =  48.13,  E 4  = 
46.12,  E5  =  44.12,  Efj  =  42.11  and  E7  =  40.11,  while  the  obtained  eigenvalues 
by  using  the  S-method  with  L  =  12  are:  Ai  =  119.40,  A2  =  50.18,  A3  =  48.19, 
A4  =  46.19,  A5  =  44.18,  A6  =  42.17,  A7  =  40.15,  A8  =  0.68, ... 

The  sensitivity  of  the  results  with  respect  to  L  is  quite  low  within  a  wide  region. 
We  have  repeated  calculations  with  values  of  L  from  L  =  10  up  to  L  =  20  and 
obtained  almost  the  same  results.  The  error  in  components  energy,  estimated  by 
corresponding  eigenvalues,  was  within  ±0.25%. 

As  a  decomposition  example,  we  omitted  the  strongest  component  and  reconstruc¬ 
ted  the  rest  of  the  signal.  The  obtained  time-frequency  representation  is  given  in  the 
upper  row  (right)  of  Figure  35.  The  last  subplot  in  Figure  35  presents  concentration 
measure  of  the  components,  in  logarithmic  scale,  which  will  be  discussed  and  used 
later. 
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Figure  35:  The  S-method  of  original  multicomponent  signal  (first  row,  left);  The 
S-method  of  reconstructed  signal  by  using  eigenvectors,  with  omitted  first 
eigenvector-component  (first  row  right);  Eigenvalues  of  the  signal’s  S-method; 
Measure  of  concentration  of  signal  components  (eigenvectors)  in  logarithmic  scale. 


46 


DRDC  Ottawa  TR  2005-251 


Figure  36:  Time-frequency  representation  of  the  eigenvectors  of  the  S-method. 
First  seven  of  them  correspond  to  the  normalized  signal  components.  In  all 
subplots  horizontal  axis  is  for  time  and  vertical  axis  is  for  frequency. 
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6  Decomposing  radar  time-varying  signals 
in  a  strong  sea-clutter 


We  now  apply  the  proposed  method  to  the  HFSWR  signals  with  strong  sea-clutter. 
Before  we  start  analysis  of  experimental  data,  we  will  consider  some  typical  situ¬ 
ations  from  the  theoretical  point  of  view. 

1 .  When  the  sea-clutter  signal  and  the  target  signal  are  separated  in  time-frequency 
plane  (occupying  unknown  and  varying  ranges  in  the  time-frequency  domain)  and 
at  the  same  time,  the  signal  energy  is  of  the  clutter  energy  order,  it  will  be  possible 
to  get  a  decomposition  such  that  the  first  component  corresponds  to  the  stronger 
clutter  and  the  second  component  corresponds  to  the  target  signal.  In  this  case 

Ui(n)  -  fciutter(ra)  =  v^ie-^Ui  (n) 

U2(n)  ~  ftargetW  =  (fl) 

where  <^>1  and  are  arbitrary  phase  constants.  Separate  time-frequency  represent¬ 
ation  of  the  clutter  and  target  signal  would  be  easy  in  this  case. 

2.  Now  consider  the  case  when  target  signal  is  much  weaker  than  the  sea  clutter 
signal  and  the  value  of  L  is  not  sufficiently  large  to  complete  the  integration  over 
clutter  and  target  signal,  according  to  (33).  Then,  first  several  components  can  come 
from  the  clutter  only  (smaller  L  will  increase  the  number  of  clutter  components). 
The  main  part  of  the  clutter  and  the  residual  clutter  parts  are  stronger  than  the  target 
signal.  Here,  we,  now,  have  to  define  a  criterion  to  select  and  resolve  the  target 
signal  component. 

3.  The  sea-clutter  signal  and  the  target  signal  are  close  or  crossing  each  other  or 
target  signal  is  highly  nonstationary  covering  wide  area.  In  this  case  both  the  clutter 
and  the  signal  could  be  separated  into  few  parts,  and  the  criterion  for  selecting  the 
target  components  should  be  defined,  as  well. 


6.1  Criterion  for  identifying  target  components 

Note  that  the  signal  components  in  HFSWR  are  single  frequency  modulated  sig¬ 
nals,  which  are  concentrated  in  the  time-frequency  plane  than  the  clutter  residual 
components.  Thus,  the  criterion  for  selecting  eigenvector(s),  being  the  target  signal 
(or  parts  of  the  target  signal),  could  be 

r  ,  s  _  (-W  +  1)  rna Xn,fc{5MU|)(w,  fc)} 

nt{P)  N/  2  y^iV/2  njyr  t  jL\ 

l^k=-N/ 2  l^n=-N/ 2  ^1V1up  \ni  K> 
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or 


CritiP ) 


(N  +  1)  ma xntk{SMUp(n,  k)} 


c\  sr'N  12  N/2 

Z  l^k=-N/2  l^n=-N/2 


SM-(n,k) 


(37) 


where  SM~p(n,  k )  denotes  negative  values  of  SMUp(n,  k ).  It  means  that  higher 
concentrated  components  (greater  maximum  in  the  numerator,  since  all  eigen¬ 
vectors  are  normalized  in  energy)  with  smaller  oscillations  (smaller  mean  absolute 
(negative)  value  in  the  denominator)  are  better  candidates  for  the  target  signal  com¬ 
ponents.  This  criterion  can  be  considered  as  a  concentration  measure,  with  infinity 
power  in  the  numerator  and  power  1  in  the  denominator  of  Crit(p).  It  can  also  be 
understood  as  the  peak-to-average  absolute  value  ratio. 


The  simplest  algorithm  is  to  find  the  value  when  Crit(p )  reaches  a  maximal  value 
(p  =  pmax)  and  assign  this  eigenvector  to  the  target  signal  component: 


ftargetH  =  Upmax(n). 

This  simple  approach  works  well  with  most  of  the  considered  signals.  Of  course, 
the  strongest  component  in  decomposition  is  the  clutter,  thus  ui  in)  is  omitted  from 
the  analysis. 

For  a  few  considered  signals,  the  target  signal  component  is  divided  over  several 
eigenvectors,  as  discussed  in  cases  2  and  3  mentioned  within  the  previous  subsec¬ 
tion.  In  order  to  deal  with  these  cases,  we  analyze  the  range  for  Cru(p)  values.  If 
the  whole  distribution  of  one  component  is  concentrated  at  one  point  in  the  time- 
frequency  plane,  i.e.,  SMUp(n,  k )  =  A8  (n  —  n0,k  —  k0),  then  Crit(p )  —  N+l.  For 
a  component  that  is  uniformly  distributed  along  a  single  line,  i.e.,  pure  linear  fre¬ 
quency  modulated  component,  Crit(p )  =  1.  For  uniformly  spread  distribution  val¬ 
ues  over  the  entire  time-frequency  plane  (obviously  not  being  a  signal  component) 
we  have  Crit (p)  =  1/(N  +  1).  Thus,  we  can  say  that  if,  for  example  Crit(i)  >  P, 
where  P  is  of  order  1  or  just  slightly  higher,  for  i  =  k,p,  then  u,  (n)  are  parts 
of  the  target  signal,  i.e., 

Crit(i)  >  P  Uj(n) is  a  target  signal  component 
^  <  P  Uj(n) is  not  a  target  signal  component 

Thus,  we  can  take  all  eigenvectors  that  satisfy  this  criterion  and  form  the  target 
signal: 

ftargetW  =  +  —  +  \/\U?(n)  (38) 

with  i  =  k,p, ...,  q  and  i  >  1,  since  the  first  component  will  always  be  the  strongest 
one  corresponding  to  the  clutter.  Note  that  if  L  in  (30)  is  not  large  enough  to 
complete  the  integration  over  the  clutter,  then  a  few  of  the  largest  components  will 
be  ones  belonging  to  the  clutter  and  i  >  2  components  should  be  omitted. 
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In  the  time-frequency  domain  we  get 

SMtaiget(n,  k )  =  g(Xk)SMk(n ,  k)  +  g(Xp)SMp(n ,  fc)  + 

...  +  g(\q)SMq(n,k)  (39) 

where  p(A)  is  a  function  of  eigenvalue.  If  we  want  to  take  all  the  components  with 
the  same  weight  g( A)  =  1,  or  if  we  want  to  keep  their  original  weights  g(  A)  =  A. 
Cases  between  these  two  are  possible,  for  example  g(X)  =  VX 

6.2  Numerical  realization 

The  numerical  realization  of  the  S-method  (30), (33)  is  very  simple,  according  to 

STFT{n ,  k)  =  DFTm{x(n  +  m)w(m)} 

and 

SM(n,k )  =  \STFT(n,k)\2  +  (40) 

L 

2  &[■ STFT(n ,  k  +  i)STFT*(n,  k  -  *)] 

2=1 

or 

SML(n,  k)  =  SML_l(n,k)  + 

2$t[STFT(n,  k  +  L)STFT*(n,  k  -  L)\ 

where  SML(n,  k )  is  SM(n,  k )  calculated  with  L  samples  in  (30)  and  SM0(n,  k )  = 
SPEC(n,  k )  =  | STFT(n,  k)  |2.  Symbol  stands  for  real  value. 

Calculation  complexity:  The  basic  step  in  the  S-method  realization  is  in  calcu¬ 
lation  of  the  STFT  by  using  the  FFT  algorithms  or  recursive  formulae  [17].  The 
calculation  of  the  STFT  in  M  time  instants  by  using  the  N  samples  FFT  requires 
an  MN  log2  N  order  of  basic  arithmetic  operations  (multiplications  and  additions). 
Additional  block  for  the  S-method  calculation,  according  to  (40),  requires  MNL 
arithmetic  operations.  For  example,  for  N  =  1024  and  L  =  10,  the  additional 
number  of  arithmetic  operations  is  of  the  same  order  as  the  basic  STFT  calculation. 
It  is  significantly  less  intensive  than  the  calculation  of  Wigner  distribution  or  any 
other  quadratic  representation  [17].  In  addition,  the  S-method  can  be  implemented 
in  hardware  [21],  which  makes  the  calculation  complexity  problem  less  important. 
Furthermore,  the  coherent  integration  time  in  the  considered  example  is  quite  long. 
The  decomposition  of  the  S-method  based  matrix  is  done  by  using  the  standard 
iterative  eigenvalue  decomposition  procedures,  with  a  given  number  of  significant 
eigenvalues  (we  used  36  significant  eigenvalues).  We  checked  the  calculation  time 
in  MATLAB  algorithms.  The  eigenvalue  decomposition  increased  this  time  for 
29%  of  the  time  required  for  the  calculation  of  the  STFT  and  the  S-method  with 
L  =  10,  N  =  1024  and  M  =  512. 
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7  Data  analysis  and  Results 


The  signals  considered  here  are  experimental  data  of  aircraft  returns,  as  used  in  [23], 
The  aircraft  is  a  King-Air  200  performing  maneuvers,  tracked  by  a  HFSWR,  with 
a  10-element  linear  receiving  antenna  array.  The  radar  carrier  frequency  is  5.672 
MHz  and  the  pulse  repetition  frequency  is  9.17762  Hz.  Each  trial  corresponds  to  a 
block  of  256  pulses.  Therefore  the  CIT  (coherent  integration  time)  of  each  signal  is 
27.89  sec.  The  King- Air  performed  two  figure-of-eight  maneuvers,  Figure  37.  Each 
figure-of-eight  maneuver  consisted  of  two  circles  with  an  approximate  diameter  of 
10km.  The  first  figure-of  eight  maneuver  was  performed  at  200  ft  (61  m),  while 
the  second  figure-of-eight  maneuver  was  performed  at  500  ft  (152  m).  The  location 
of  the  King-Air,  when  each  signal  was  collected,  is  marked  by  a  square,  Figure 
37.  Each  signal  represents  a  different  scenario  that  could  arise  when  tracking  a 
maneuvering  aircraft. 

In  the  HFSWR,  beside  the  target,  the  signal  contains  clutter  that  is  primarily  due 
to  scattering  from  the  surface  of  the  sea.  The  received  signal  of  a  pulse  Doppler 
radar,  taken  at  a  particular  range,  is  viewed  in  the  Fourier  transform  or  Doppler 
domain.  Two  sharp  peaks,  Bragg  components,  can  be  observed.  These  peaks  indic¬ 
ate  that  the  dominant  form  of  sea-clutter  is  due  to  scattering  from  sea  waves,  with 
wavelengths  being  half  of  the  wavelength  of  the  radar  carrier  frequency,  travelling 
radially  towards  and  away  from  the  radar  [7],  [1].  The  Bragg  components  of  the 
clutter  are  called  1st  order  clutter.  In  addition  to  the  Bragg  components,  the  Doppler 
spectrum  has  a  continuum  called  the  2nd  order  clutter. 

7.1  Calculation  procedure: 

Denote  signal  length  by  N.  In  all  considered  cases  we  use  N  =  256. 

Step  1 .  Calculate  the  STFT  (with  rectangular  window)  of  the  zero-padded  signal, 
oversampled  by  factor  2. 

Step  2.  Calculate  the  S-method  of  the  signal  according  to  (40)  for  a  given  L,  for 
example  L  =  10.  Size  of  obtained  S-method  is  2N  =  512  samples  in  time  and 
4 N  =  1024  samples  in  frequency.  The  central  part  of  the  S-method  in  the  frequency 
domain  is  equivalent  to  the  Wigner  distribution.  Note  that  the  frequency  range  of 
the  S-method  is  like  in  the  Fourier  transform  (— /s/2,  fs/ 2)  and  (— /s/4,  fs/ 4)  for 
the  Wigner  distribution,  where  fs  is  the  sampling  frequency.  Only  even  samples  are 
included  in  further  analysis  in  order  to  avoid  non-integer  indices  in  (35). 

Step  3.  Calculate  matrix  R  according  to  (35).  Since  the  signal  is  zero-padded,  use 
only  even  rows  in  order  to  avoid  non-integer  indices.  Since  only  central  part  of 
S-method,  and  only  even  samples  are  used,  order  of  the  matrix  R  is  N  =  256 
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Figure  37:  Path  of  the  King-Air  200  as  a  function  of  range  (in  km)  and  azimuth  (in 
degrees) 


Step  4.  Decompose  R  into  eigenvectors  and  eigenvalues.  Note  that  eigenvectors 
has  the  same  length  as  the  original  signal. 

Step  5.  Calculate  the  S-method  (40)  of  the  eigenvectors  by  using  the  STFT  with  a 
window,  for  example  square  root  Hanning  window  of  256  samples. 

Step  6.  Calculate  measure  (37)  and  assign  the  eigenvector  corresponding  to  the 
maximal  measure  to  the  target  signal  if  C~it(jpmax)  >  2.  If  C~it  (pmax)  <  2  then 
there  is  no  target  signal  detected  (target  signal  is  too  close  to  the  clutter).  In  this 
case,  repeat  steps  1.  to  6.  with  smaller  L,  for  example  L  —  4.  If  C~t(pm ax)  > 
2  check  if  there  is  any  other  highly  concentrated  eigenvector,  for  example  with 
C~it(p  7^  Pmax)  >  3.  If  there  is  such  a  highly  concentrated  eigenvector,  it  should 
also  be  a  part  of  the  target  signal.  Include  it  according  to  (38). 

Step  7.  Show  the  time-frequency  representation  of  the  resulting  target  signal  (44) 
and  calculate  its  high  resolution  version  (45),  with  excluded  isolated  points.  High 
resolution  images  are  obtained  by  using  the  high  resolution  version  of  the  S-method 
presented  in  the  Appendix.  After  a  high  resolution  distribution  is  calculated,  the 
pattern  recognition  algorithm  is  used  to  eliminate  instants  when  isolated  points  are 
produced  (random  maxima  at  the  instants  when  there  is  no  target  signal  compon¬ 
ent).  Only  the  regions  where  more  than  10  connected  points  of  the  high  resolution 
S-method  exist,  are  kept  for  final  high  resolution  presentation. 

Step  8.  Take  the  next  signal  and  go  to  Step  1. 
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By  using  this  procedure  we  analyze  all  positions  for  the  described  experiment.  The 
target  signal  is  detected  and  separated  in  all  cases.  In  four  cases,  the  target  signal 
was  too  close  to  the  clutter,  and  the  target  was  not  detected  when  L  =  10.  The 
procedure  was  repeated  with  L  =  4,  according  to  the  description  in  step  6,  and  the 
target  was  then  detected. 

Now,  the  proposed  procedure  is  used  to  decompose  the  signal.  The  following  typ¬ 
ical  cases  are  presented  here. 

•  A  stationary  target  signal  (aircraft  is  moving  with  a  constant  velocity),  far  from 
the  clutter,  Figure  38.  The  time-frequency  representation  of  the  corresponding 
eigenvectors  is  shown  in  Figure  39. 

•  A  nonstationary  target  signals  far  from  the  clutter,  Figures  40,  42  along  with 
time-frequency  representation  of  the  eigenvectors,  Figure  41,  for  the  signal 
shown  in  Figure  40. 

•  A  highly  nonstationary  target  signal  covering  wide  frequency  range,  Figure  43. 

•  A  nonstationary  target  signal  very  close  to  the  clutter,  Figure  44. 

•  A  nonstationary  target  signal  intersecting  the  clutter,  Figure  45. 

For  all  the  signals,  the  time-frequency  representation  (S-method)  of  the  original  sig¬ 
nal  including  the  clutter,  is  given  (upper  right  subplot  in  Figures  38,40,42,43,44,45). 
Since  the  clutter  is  extremely  strong,  in  order  to  get  signal  component  visible,  the 
limiter  is  used  in  time-frequency  representation.  The  upper  left  subplot  represents 
concentration  measure  for  eigenvectors  in  all  considered  cases.  Lower  subplots 
present  the  S-method  of  the  detected  target  signal  component  (right)  and  its  highly 
concentrated  version  (left).  The  frequency  axis  in  the  highly  concentrated  repres¬ 
entations  is  scaled  to  represent  target  instantaneous  radial  velocity  while  dashed 
line  represents  target  radial  velocity  obtained  by  standard  FFT  technique.  Figures 
39  and  41  present  time-frequency  representations  of  the  eigenvectors  used  for  the 
concentration  measure  and  target  signal  component  detection.  Plots  for  the  whole 
trial  are  presented  in  appendix  2. 

For  comparison  with  the  FFT  method,  three  cases  are  presented:  non-accelerating 
target  far  from  the  Bragg’s  lines,  accelerating  target  far  from  the  Bragg’s  lines  and 
the  target  very  close  to  the  Bragg’s  lines,  Figure  46.  As  expected,  in  the  case  of 
constant  velocity  (non-accelerating  target)  the  FFT  method  produces  a  clear  result 
with  a  peak  corresponding  to  the  Doppler  frequency.  It  is  in  accordance  and  in  good 
agreement  with  the  time-frequency  result  presented  in  Figure  38.  However,  in  the 
second  case  the  FFT  method  is  smeared,  and  the  constant  frequency  estimated  based 
on  the  FFT  does  not  correspond  to  the  real  event  of  fast  varying  target  velocity, 
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Figure  38:  Signal  3  (reference  to  Figure  1.):  Time-frequency  representation  (top 
right),  Concentration  of  the  eigenvectors  time-frequency  representations  (top  left), 
Time-frequency  representation  of  the  detected  target  signal  (bottom  right),  Highly 
concentrated  time-frequency  representation  of  the  detected  target  signal  (bottom 
left)  and  FFT  estimation  of  the  target  velocity  (dashed  line) 


Figure  43.  In  the  third  case,  there  is  a  slightly  smeared  peak,  but  connected  to  the 
clutter  spectrum.  Using  the  FFT  method  it  would  be  difficult  to  conclude  that  a 
target  exists.  The  time-frequency  based  approach  clearly  indicates  that  we  have  a 
target  signal,  including  its  separation  by  the  proposed  method,  Figure  44.  Note  that 
the  velocities  obtained  by  using  the  FFT  approach  are  presented  for  all  considered 
cases,  as  mentioned  earlier. 
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Figure  39:  Time  frequency  representations  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  signal  detection  for  signal  presented  in  Figure 
5.  In  all  subplots  horizontal  axis  is  for  time  and  vertical  axis  is  for  frequency. 
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Figure  40:  Signal  31:  Time-frequency  representation  (top  right),  Concentration  of 
the  eigenvectors  time-frequency  representations  (top  left),  Time-frequency 
representation  of  the  detected  target  signal  (bottom  right),  Highly  concentrated 
time-frequency  representation  of  the  detected  target  signal  (bottom  left)  and  FFT 
estimation  of  the  target  velocity  (dashed  line) 
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Figure  41:  Time  frequency  representations  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  signal  detection  for  signal  presented  in  Figure 
7.  In  all  subplots  horizontal  axis  is  for  time  and  vertical  axis  is  for  frequency. 
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Figure  42:  Signal  43:  Time-frequency  representation  (top  right),  Concentration  of 
the  eigenvectors  time-frequency  representations  (top  left),  Time-frequency 
representation  of  the  detected  target  signal  (bottom  right),  Highly  concentrated 
time-frequency  representation  of  the  detected  target  signal  (bottom  left)  and  FFT 
estimation  of  the  target  velocity  (dashed  line) 
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Figure  43:  Signal  53:  Time-frequency  representation  (top  right),  Concentration  of 
the  eigenvectors  time-frequency  representations  (top  left),  Time-frequency 
representation  of  the  detected  target  signal  (bottom  right),  Highly  concentrated 
time-frequency  representation  of  the  detected  target  signal  (bottom  left)  and  FFT 
estimation  of  the  target  velocity  (dashed  line) 
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Figure  44:  Signal  7:  Time-frequency  representation  (top  right),  Concentration  of 
the  eigenvectors  time-frequency  representations  (top  left),  Time-frequency 
representation  of  the  detected  target  signal  (bottom  right),  Highly  concentrated 
time-frequency  representation  of  the  detected  target  signal  (bottom  left)  and  FFT 
estimation  of  the  target  velocity  (dashed  line) 


60 


DRDC  Ottawa  TR  2005-251 


15 


Components  concentration 


S-  method 


<j> 

c 

o 

) 

0  10  20  30 


eigenvector  No 


time 


Targetsignal  HC  S-method  of  target  signal 


Figure  45:  Signal  19:  Time-frequency  representation  (top  right),  Concentration  of 
the  eigenvectors  time-frequency  representations  (top  left),  Time-frequency 
representation  of  the  detected  target  signal  (bottom  right),  Highly  concentrated 
time-frequency  representation  of  the  detected  target  signal  (bottom  left)  and  FFT 
estimation  of  the  target  velocity  (dashed  line) 
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Figure  46:  Fourier  transform  (aboslute  value)  of  the  three  typical  signals: 
nonaccelerating  target  far  from  Bragg’s  lines,  accelerating  target  far  from  Bragg’s 
lines  and  target  very  close  to  Bragg’s  lines.  Triangle  arrows  show  target  signal. 
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8  Noise  Analysis  in  the  Algorithm  for 
Signal  Decomposition 


Let  us  consider  four  component  signal,  where  each  component  is  a  Gaussian  chirp, 
contaminated  with  a  white  Gaussian  complex  noise  e(n)\ 

4 

x(n)  =  ^Akej^n-dk)+jak{-^^e-(J^^  +  e(n) 
k= 1 

where  d4  =  d2  =  —64,  d3  =  0,  d4  =  64,  ,  uj1  =  u>4  =  0,  u>2  =  —  oj 3  = 
a i  —  a:i  =  2^:,  a2  =  a4  =  —  and  —128  <  n  <  127.  Note  that  the  signal 
components  are  separated  in  the  time-frequency  plane. 

In  the  noiseless  case  we  will  assume  that  e(n)  =  0  and  A\  =  1.3,  ^42  =  1.2, 
^3  =  1.1  and  A4  =  1.0.  Rectangular  window  of  64  samples  length  is  used  and 
L  =  36  is  chosen  in  order  to  satisfy  conditions  form  Proposition  1 .  The  results 
are  presented  in  Table  1  and  Figure  47.  The  energy  of  each  signal  component  and 
corresponding  eigenvalue  is  presented  in  Table  1.  Note  that  eigenvalues  correspond 
well  to  the  components  energies.  Figure  47  presents  spectrogram  of  the  original 
signal,  TFRs  of  the  first  four  eigenvectors  and  TFR  of  the  reconstructed  signal. 


Component 

1 

2 

3 

4 

Energy 

67.8 

57.7 

48.5 

40.1 

Eigenvalue 

67.5 

57.5 

48.3 

39.9 

Table  1:  Component  energies  and  corresponding  eigenvalues 


The  decomposition  algorithm  is  applied  to  the  considered  signal  with  =  ^42  = 
A3  =  A4  =  1  for  signal  to  noise  ratio  0  dB  (Figure  48)  and  —4  dB  (Figure  49). 
Hanning  window,  128  samples  length  and  L  =  16  is  used  in  booth  cases.  In  both 
cases  all  signal  components  are  separated,  and  the  reconstructed  signal  is  obtained 
without  noise  in  the  parts  of  the  time-frequency  plane  where  there  are  no  signal 
components.  In  noisy  cases,  equal  components  energies  can  be  used  because  high 
noise  introduce  different  eigenvalues  in  the  decomposition  process,  avoiding  pos¬ 
sible  ambiguity. 


8.1  Noise  analysis 

If  the  analyzed  signal  is  corrupted  with  noise,  then  the  assumption  that  only  M 
eigenvalues  are  different  from  zero  is  not  valid.  On  the  other  hand,  the  noise  com¬ 
ponents  are  distributed  over  all  N  eigenvectors. 


62 


DRDC  Ottawa  TR  2005-251 


The  described  decomposition  algorithm  is  analyzed  for  various  SNR  and  with  dif¬ 
ferent  algorithm  parameters  (L  in  the  S-method  calculation  and  window  length  in 
the  STFT  calculation).  The  results  are  shown  in  Tables  2  and  3.  For  each  combin¬ 
ation  of  SNR,  L  and  window  length  h,  the  algorithm  is  repeated  100  times,  giving 
total  of  400  components  for  detection.  The  percentage  of  missed  components  is 
shown  in  tables  and  used  as  measure  of  algorithm  robustness. 


SNR 

L  =  8 

L  =  16 

L  =  32 

L  =  64 

OdB 

3.50% 

0.00% 

0.00% 

0.00% 

-2dB 

6.50% 

0.00% 

0.50% 

0.00% 

-4dB 

10.00% 

3.25% 

2.50% 

3.50% 

-6dB 

20.25% 

14.75% 

11.50% 

7.50% 

-8dB 

35.50% 

29.50% 

22.75% 

21.00% 

Table  2:  Sensitivity  of  the  proposed  algorithm  to  the  choice  of  L  for  various  SNR. 
Window  length  is  128  samples. 


SNR 

h  =  160 

h  =  128 

h  =  64 

OdB 

0.25% 

0.00% 

0.50% 

-2dB 

0.50% 

0.50% 

0.75% 

-4dB 

3.50% 

2.50% 

4.50% 

-6dB 

16.25% 

11.50% 

14.00% 

-8dB 

28.00% 

22.75% 

33.25% 

Table  3:  Sensitivity  of  the  proposed  algorithm  to  the  choice  of  window  length  h  for 
various  SNR  with  L  =  64 

It  is  shown  that  decomposition  algorithm  is  very  robust  with  respect  to  the  para¬ 
meter  values.  The  heuristic  analysis  gives  some  estimations  of  parameter  values. 
Namely  if  we  want  to  obtain  Wigner  distribution  then  the  parameter  L  should  be 
equal  to  the  half  of  the  component  frequency  support  (in  discrete  domain).  Time 
window  in  STFT  calculation  should  be  long  enough  so  the  whole  component  is 
covered  by  window.  On  the  other  hand  very  large  window  can  combine  two  com¬ 
ponents  in  one  eigenvector. 
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Spectrogram  Eigenvector  No  1  Eigenvector  No  2 


time  time  time 


Figure  47:  Spectrogram,  reconstructed  signal  and  TFRs  of  the  first  four 
eigenvectors  -  noiseless  case 
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Figure  48:  Spectrogram  of  noisy  signal,  eigenvectors  TFRs  and  TFR  of  the 
reconstructed  signal  for  OdB  SNR 
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Figure  49:  Spectrogram  of  noisy  signal,  eigenvectors  TFRs  and  TFR  of  the 
reconstructed  signal  for  —4 dB  SNR. 
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Figure  50:  Decompostion  of  a  real  HF  radar  signal  in  a  strong  sea-clutter 
(realization  1) 
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Figure  51:  Decompostion  of  a  real  HF  radar  signal  in  a  strong  sea-clutter 
(realization  2) 
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Figure  52:  Decompostion  of  a  real  HF  radar  signal  in  a  strong  sea-clutter 
(realization  3) 
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9  Conclusion 


The  analysis  and  processing  of  radar  signals  are  usually  performed  in  the  time  or 
frequency  domain.  These  methods  can  be  inefficient  in  the  cases  when  significant 
time  variations  of  frequency  content  are  present,  like  as  in  the  High-Frequency 
Surface-Wave  Radar  (HFSWR)  radar  signals  considered  in  this  report. 

This  report  examines  the  possibility  of  improving  the  detection  of  accelerating 
targets  using  High-Frequency  Surface- Wave  Radar  (HFSWR)  in  the  presence  of 
sea-clutter.  Because  we  are  dealing  with  high  speeds  and  accelerations,  the  time- 
frequency  analysis  has  the  capability  of  improving  detection,  since  its  resolution 
abilities  are  not  reduced  by  a  spread  in  the  Doppler  signature. 

In  this  report,  we  proposed  a  new  time-frequency  based  signal  decomposition 
method.  This  approach  is  based  upon  the  eigenvalue  decomposition  method,  com¬ 
bined  with  the  S-method  that  produces  a  sum  of  the  Wigner  distribution  of  indi¬ 
vidual  signal  components.  The  proposed  decomposition  method  is  theoretically 
derived.  The  efficiency  and  accuracy  of  the  proposed  decomposition  method  is 
demonstrated  on  simulated  examples.  Then,  it  is  used  in  improving  the  analysis 
of  signals  obtained  from  a  maneuvering  low-altitude  aircraft  using  HFSWR  in  the 
presence  of  sea-clutter.  Since  we  are  dealing  with  high  speeds  and  accelerations 
of  the  aircraft,  the  time-frequency  method  improves  the  detection,  because  its  res¬ 
olution  abilities  are  not  reduced  by  a  spread  in  the  Doppler  signature.  The  study 
performed  by  using  the  new  method  demonstrates  that  the  use  of  time-frequency 
signal  decomposition  can  improve  target  velocity  estimation  and  detection  perform¬ 
ance  of  the  HFSWR.  The  method  provides  true  time-varying  Doppler  shift  within 
the  considered  time,  which  is  not  possible  with  the  Fourier  transform  method  that 
provides  only  its  average  value.  In  addition,  the  proposed  method  successfully  de¬ 
tects  the  target  signal  in  all  cases  of  the  presented  real  experiment,  which  is  not  the 
case  when  the  detection  is  performed  by  the  Fourier  transform  method. 

The  noise  analysis  is  performed  in  the  algorithm  for  signal  decomposition  and  it 
is  shown  that  the  choice  of  algorithm  parameters  does  not  have  a  great  influence 
on  the  decomposition  process.  The  method  presented  here  is  not  restricted  to  this 
application,  but  it  can  be  applied  also  in  various  other  settings  of  nonstationary  sig¬ 
nal  analysis  and  filtering,  for  example,  moving  target  identification  in  Inverse  Syn¬ 
thetic  Aperture  Radar  (SAR)  images  and  micro-Doppler  separation  from  Inverse 
Synthetic  Aperture  Radar  (ISAR)  images,  etc. 

Although  this  new  approach  demonstrates  the  improved  detection  performance  of 
an  accelerated  air  target  using  HFSWR  in  sea-clutter  environment,  a  more  detailed 
analysis  will  be  performed  in  a  future  study  in  relation  to  other  clutter  environments 
such  as  ionosphere,  meteors,  and  lighting  environments.  Furthermore,  the  s-method 
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based  detectors  will  be  developed  for  real-time  operational  scenarios  and  will  then 
be  compared  with  existing  Fourier  Transform  based  detectors.  Once  the  new  ap¬ 
proach  outperforms  the  currently  available  detector  methods,  numerically  efficient 
algorithms  will  be  developed  for  real-time  operations. 
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10  Appendix  1:  High-resolution  S-method 
with  relation  to  the  Capon’s  method 


As  shown  in  [16]  all  spectral  estimators  can  be  considered  as  either  being  smoothed 
or  being  cross-calculated.  The  first  family  leads  to  the  variations  of  smoothed  spec¬ 
trograms,  while  the  other  leads  to  the  variations  of  the  estimators  defined  in  [17]. 
Here,  we  will  show  that  the  Capon’s  method  can  be  considered  as  a  version  of  the 
smoothed  spectrograms  (spectrograms  of  eigenvectors  for  close  signal  compon¬ 
ents  and  spectrograms  of  signal  components  for  separated  components).  A  form 
of  highly  concentrated  distributions  based  on  the  S-method  may  be  defined  in  an 
analog  manner,  using  the  analysis  presented  in  this  paper. 


The  Capon’s  filtering  method  applied  to  N  +  1  samples  of  a  signal  f(n),  being  a 
sum  of  complex  sinusoids  exp^'u^n),  denoted  in  a  vector  form  as  f  (n)results  in  the 
distribution 


Scapij^i 


1 

a*(A:)Rfr1a(fc)  ’ 


(41) 


where  *  denotes  Hermitian  transpose  and  a (k)  =  [1  ej‘27rk/(N+1)  ...e?27rfe]T,  with 


Rf  =  £{f(n)f*(n)}. 


(42) 


In  practice,  the  autocorrelation  matrix  Rf  is  estimated  by: 


K 

Rf  (n,K)  =  ±  ^f(n  +  p)f>  +  p)  +  pI=£QQ*+pI,  (43) 

p= i 


where  I  is  the  identity  matrix  used  for  regularization  and  Q  is  the  matrix  whose 
columns  are  signal  vectors  f (n  +  p),p—  1, ...,  K. 


The  Capon’s  form  can  be  written  by  using  eigenvalue  decomposition  of  the  auto¬ 
correlation  matrix  Rf  (n,  K)  as 

Rf  (n,  K)  =  -^QQ*  +  pl  =  =  ^VAV*  +  pi, 

By  using  this  decomposition,  we  can  write 

i-  D1  +  tVv„v; 

i  ^ p 

p=  1  1  J 

where  V;,  are  eigenvectors  and  \p  are  eigenvalues  of  Rf(n,  K )  The  Capon’s  form 
then  reads 

S C apK  ijl  >  > 

N+l-Y,K%rrSPECVr(n,k) 

P=  1 


Rf“1(^^)  =  - 

P 
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where  SPECyp(u,  t)  are  the  spectrograms  of  the  eigenvectors  and  p  — >  0.  For  sep¬ 
arated  components  the  eigenvectors  correspond  to  the  signal  components,  summed 
with  the  same  weights. 

In  our  case,  the  S-method  form  that  corresponds  to  the  Capon’s  form  would  lead  to 
the  factor  of 

K 

SMnor(n,  k )  =  ^  SMUp(n,  k ) 

P=  1 
K 

instead  of  Y,  SPECvp(n ,  k)  in  the  Capon’s  form,  resulting  in  a  distribution: 

p= i 

SMCK(n,k )  =  - ^ -  (44) 

iV  +  l-£ppSMu>,,fc) 

p=i 

A 

_  K  ’ 

1  -  E  cpSMUp  (n,  k) 

P=  1 

where  A  =  p/(iV  +  1)  and  cp  =  pp/ (N  +  1). 

In  the  realization  we  used  a  very  simple  normalized  version  of  the  highly  concen¬ 
trated  S-method: 

SMCK(n,  k )  = - - - 1 - .  (45) 

■A  1.01maxfe{SM„  (n,k)} 

n — 1  V 
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1 1  Appendix  2:  Plots  for  the  whole  trial 
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Figure  53:  Signal  1-  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  54:  Signal  1  -  The  frequency  representation  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is  for 
time  and  vertical  axis  is  for  frequency 


DRDC  Ottawa  TR  2005-251 


75 


S-method  S-method  of  target  signal 


Components  concentration 


0  10  20  30 


Target  signal  HC 


time 


eigenvector  No 


Figure  55:  Signal  2  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  56:  Signal  2  -  The  frequency  representation  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is  for 
time  and  vertical  axis  is  for  frequency 
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Figure  57:  Signal  3  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  58:  Signal  3  -  The  frequency  representation  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is  for 
time  and  vertical  axis  is  for  frequency 
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Figure  59:  Signal  4  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  60:  Signal  4  -  The  frequency  representation  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is  for 
time  and  vertical  axis  is  for  frequency 
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Figure  61:  Signal  5  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  62:  Signal  5  -  The  frequency  representation  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is  for 
time  and  vertical  axis  is  for  frequency 
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Figure  63:  Signal  6  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  64:  Signal  6  -  The  frequency  representation  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is  for 
time  and  vertical  axis  is  for  frequency 
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Figure  65:  Signal  7  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  66:  Signal  7  -  The  frequency  representation  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is  for 
time  and  vertical  axis  is  for  frequency 
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Figure  67:  Signal  8  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  68:  Signal  8  -  The  frequency  representation  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is  for 
time  and  vertical  axis  is  for  frequency 
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Figure  69:  Signal  9  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  70:  Signal  9  -  The  frequency  representation  of  the  eigenvectors  used  for  the 
concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is  for 
time  and  vertical  axis  is  for  frequency 
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Figure  71:  Signal  10  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  72:  Signal  10  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  73:  Signal  1 1  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  74:  Signal  11  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  75:  Signal  12  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  76:  Signal  12  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 


DRDC  Ottawa  TR  2005-251 


97 


S-method  S-method  of  target  signal 


Components  concentration 


Target  signal  HC 


■K 

i— 

■E 

O 

C 

o 

o 


time 


Figure  77:  Signal  13  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  78:  Signal  13  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  79:  Signal  14  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  80:  Signal  14  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  81:  Signal  15  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  82:  Signal  15  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  83:  Signal  16  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  84:  Signal  16  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  85:  Signal  1 7  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  86:  Signal  17  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  87:  Signal  18  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  88:  Signal  18  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  89:  Signal  19  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  90:  Signal  19  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  91:  Signal  20  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  92:  Signal  20  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  93:  Signal  21  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  94:  Signal  21  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  95:  Signal  22  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  96:  Signal  22  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 


DRDC  Ottawa  TR  2005-251 


117 


time 


S-method  of  target  signal 


>■ 

o 

c 

rv 

■i= 


time 


c 

g 

l— 

■E 

O 

c 

o 

o 


Components  concentration 
30 

20 

10 

0 

0  10  20  30 


< 

■f 

Target  signal  HC 


eigenvector  No 


Figure  97:  Signal  23  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  98:  Signal  23  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  99:  Signal  24  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  100:  Signal  24  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  101:  Signal  25  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  102:  Signal  25  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  103:  Signal  26  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  104:  Signal  26  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  105:  Signal  27  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  106:  Signal  27  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  107:  Signal  28  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  108:  Signal  28  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  109:  Signal  29  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  110:  Signal  29  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  111:  Signal  30  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  112:  Signal  30  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  113:  Signal  31  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  114:  Signal  31  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  115:  Signal  32  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  116:  Signal  32  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  117:  Signal  33  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  118:  Signal  33  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  119:  Signal  34  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  120:  Signal  34  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  121:  Signal  35  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  122:  Signal  35  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  123:  Signal  36  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  124:  Signal  36  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  125:  Signal  37  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  126:  Signal  37  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  127:  Signal  38  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  128:  Signal  38  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  129:  Signal  39  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  130:  Signal  39  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  131:  Signal  40  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  132:  Signal  40  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  133:  Signal  41  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  134:  Signal  41  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  135:  Signal  42  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  136:  Signal  42  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  137:  Signal  43  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  138:  Signal  43  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  139:  Signal  44  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  140:  Signal  44  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  141:  Signal  45  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  142:  Signal  45  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  143:  Signal  46  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 


164 


DRDC  Ottawa  TR  2005-251 


Figure  144:  Signal  46  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  145:  Signal  47  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  146:  Signal  47  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  147:  Signal  48  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  148:  Signal  48  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  149:  Signal  49  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  150:  Signal  49  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  151:  Signal  50  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 


172 


DRDC  Ottawa  TR  2005-251 


Figure  152:  Signal  50  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  153:  Signal  51  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  154:  Signal  51  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  155:  Signal  52  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  156:  Signal  52  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  157:  Signal  53  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  158:  Signal  53  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  159:  Signal  54  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  160:  Signal  54  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  161:  Signal  55  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  162:  Signal  55  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  163:  Signal  56  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  164:  Signal  56  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  165:  Signal  57  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  166:  Signal  57  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  167:  Signal  58  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  168:  Signal  58  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  169:  Signal  59  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  170:  Signal  60  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  171:  Signal  60  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  172:  Signal  60  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  173:  Signal  61  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  174:  Signal  61  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  175:  Signal  62  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  176:  Signal  62  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  177:  Signal  63  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  178:  Signal  63  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  179:  Signal  64  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  180:  Signal  64  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  181:  Signal  65  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  182:  Signal  66  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  183:  Signal  66  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  184:  Signal  67  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  185:  Signal  67  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  186:  Signal  67  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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Figure  187:  Signal  68  -  Concentration  of  the  eigenvectors  time-frequency 
representations  (bottom  left),  Time-frequency  representation  of  the  detected  target 
signal  (top  right),  Highly  concentrated  time-frequency  representation  of  the 
detected  target  signal  (bottom  right) 
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Figure  188:  Signal  68  -  The  frequency  representation  of  the  eigenvectors  used  for 
the  concentration  calculation  and  target  detection.  In  all  subplots  horizontal  axis  is 
for  time  and  vertical  axis  is  for  frequency 
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